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We use the method of QCD sum rules to investigate the isospin symmetry breaking of K and 
K* mesons. The electromagnetic effect, difference between up and down current-quark masses and 
\^ • difference between up and down quark condensates are important. We perform sum rule analyses 

' of their masses and decay constant differences, which are consistent with experimental values. Our 

\ results yield A/x = f^a - = 1.5 MeV. 
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' QCD has an approximate flavor symmetry which is determined by the pattern of the quark masses. Isospin 

^sO ', symmetry in particular holds to a high accuracy. This is because the scale is set by (to,; — md)/A^, where m's 

■ are current quark masses, while A^^ is the chiral symmetry breaking scale around 1 GeV. Because of this small 

' hadronic isospin violations, the electromagnetic effect becomes important in order to understand the isospin symmetry 

, breaking JJJa, Q . There are many papers suggesting that the electromagnetic effect is dominant in the mass splitting 

' of pions 0, S • 

Qh 

I ' Therefore, to study the isospin symmetry breaking, it is necessary to consider both the hadronic isospin violations 

Oh! and the electromagnetic effect. In this paper, we study the isospin symmetry breaking of K and K* [J^ — 0+) 

^ mesons in the QCD sum rule. This work is an extension of the previous one for the tt and p mesons @. 

^ I One can calculate the hadronic effect due to the different up and down current-quark masses and condensates. While 
•'H , for the electromagnetic effect, we follow the procedure in Ref. j2|]. They constructed a gauge invariant electromagnetic 
' two-point function for the heavy-light quark systems. 

H ; . . . 

^ This paper is organized as follows. In section 2, we derive the QCD sum rules for the K and K* mesons. In section 
3, we discuss our numerical results of their masses, decay constants and differences. We find that they are consistent 
with the experimental values. Section 4 is a summary. 



I. INTRODUCTION 
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II. QCD SUM RULES FOR K AND K* MESONS 

For the past decades QCD sum rule has proven to be a powerful and successful non-perturbative method [sL IT(tI|. 
In sum rule analyses, we consider two-point correlation functions: 



(1) 



where for K meson 



and for the vector K* meson 



»7^^' = 917/^7592 . 



(2) 



(3) 



Here these currents may couple to particles K and K* through 

(0|r?f' )(0)|if*(p)) = fK'mK^4 



K* 



Here is the four momentum carried by the initial meson, and fK- are the decay constants of K and K* 
respectively, niK* is the mass of if*, and is the polarization vector of K* . In the OPE, H(q^) can be divided into 
two parts: the hadronic part and the contributions from the electromagnetic effects. The hadronic part for K and 
K* have been calculated in the original work of the QCD sum rule 0,U 

For the charge neutral current, like and K*^ , we can change the gluons in QCD to the photons up to the order 
of Q!e (= e^/47r), and easily calculate electromagnetic contributions. For the charged current, like and K*^^ the 
calculation of electromagnetic contributions is slightly more complicated. If we simplychange gluons to photons, the 
result is not gauge invariant. To solve this problem, we follow the procedure in Ref. Expanding to order cte, the 
currents become (Fig. 



917^7592 - eer qi{q - h) 7^75 



(fci - ^2)1/ 
(fci - fc2)2 



92(^2) : 



9i7m'?2 - eer - fci) 7^ 



(fci 



(fci - 

where the normalized total charge of the meson is defined by er = Sqi 



A^{ki - k2) (J2(fc2) , 

- e„2, and takes ±1 for K"^ and K*^ . 



(4) 
(5) 






FIG. 1: The gauge invariant current up to order 



We have performed the OPE calculation up to dimension six, which contains the four-quark condensates. The 
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results are 



^Kf'(KO) - ^TTV-'^^ — 2 o 2~ 



3 1 _ 
— — (TOrf + TO^)^ In — -{md{ss) +ms{dd)) (6) 

1 ,Q!., 2\ 3527r ,T,, 447r nJj\2 . \2\ 

-G^) - -—^as{dd){ss) - —^aeedesiiddy + (ssy) , 



-, , „-(to,i+TOs) In — 5 j{mu{ss) + ms{uu)) (7) 

1 ,a.s^2\ 3527r 447r 2 , /- \2\ 

Lr ) - — — 7rQ!s(uu)(ss) ~ -—-T:aeeues{{uu) + (ss) ) , 



ii/^.of !(-*o\ — — - — ^(-L H Jin — — m — — 



3 -q^ 1 _ 

' 099 ("^rf ~ In — 5- + -j{md{ss) + ms{dd)) (8) 

'^^aeedes{{ddf + (ss)^) , 



TT 1 /I , 1 1 

47r^ TT /x^ 47r-^ tt /i^ 

3 -g^ 1 _ 

;{mu-ms)'^\n — 5- + — (m„(ss) + ms(uw)) (9) 



87r2g2 y^2 ^. 

1 ^ t; ^-,0 \ 327r ,_ 327r //_ \9 /_\9\ 

■ G^> + — rQ;5(uu)(ss) - — — rQ;s((ww) + (ss)'') 



12g4^7r ' 9q^ ' " ' 8lq^ 
287r 
27g' 



gaeej.(uu)(ss) + ■;^;^aee„es((uM) + (ss) ) 



In these equations, u, d, s represent u, d and s quarks respectively. The couplings Cu, Cd and are normalized by the 
unit electric charge e, and therefore, e„ = 2/3 and 6^ = 6^ = —1/3. The quantities {uu), (dd) and (ss) are dimension 
D = 3 quark condensates, and (g'^GG) is a D = 4 gluon condensate. We have assumed the vacuum dominance and 
factorization for the four quark condensates, for instance 

A° A" 16 

(OlwsyWT^TsyglO) = -{0\qq\Of , 

\a 2g 

(0|9a^.75YOT^^''75Y9|0) - y(0|gg|0)^ 



The difference Hk" ~ determines the isospin symmetry breaking of K meson, while the difference 11^. — Ilj^t± 
determines the isospin symmetry breaking of K* meson. If we consider that the difference between the up and down 
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quark condensates is small and introduce the average condensate (qq) = {{uu) + {dd))/2, we find 



1 , 2 2\Q^ei —q^ 3 , r, \/ 

— ^(2e„es - - ej — In ^ (mu +md + 2ms)[m,u - md) In — ^ 



+ \{mu - md){ss) + -^ms{{uu) - {dd)) 

3527]- _ 4471" 

+ -^Y^o;s{ss)i{uu) - (dd)) + ^y^ae(e„ - ed)es{{qqf + {ssf) , 

1 /o 2 2x"e-, -9^ 3 
— r(2e„es - - eJ — In — (niu +md- 2ms)(m„ - md) In — r- 



-\{mu - md){ss) - \ins{{uu) - (dd)) - ^^as{ss){{uu) ~ (dd)) 

9q" 



+ ^^Ois{{uuf - (ddf) - ^aeeT{uu){ss) - - ed)es{{qqf + {ssf) . 

There are three non-perturbative effects 

1. The difference due to the masses of up and down quarks. 

2. The difference between the up and down quark condensates. 

3. The electromagnetic part containing four-quark condensates which are of the first order of ag. 

The difference between the up and down quark condensates has been evaluated previously. We define A to be 

A ^ 7^ - 1 10 

\uu) 

For instance, Gasser and Leutwyler obtained A « —0.0074 while in Ref Hatsuda, Hogaasen and Prakash 
found -0.0078 < A < -0.0067. In the QCD sum rule, Chernyak and Zhitnitsky obtained A « -0.009 [3. Here we 
will use the value A w -0.0074. 

If we choose <f ~ m|-, the above three effects are in the same order of magnitude. This is different from the tt and 
p mesons, where only the electromagnetic part dominates Q. 

Within the approximation of the narrow resonance with a continuum above threshold value sq, after the Borel 
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transformation, we obtain the final QCD sum rules 

= + — + e<je, — )Mi(l-e"^)- A(™<i + m,)2lnAfB(l-e"^) 

47r^ TT TT 47r^ 

1 _ 1 /T 

+ m,{dd)) + ^^(-IG^) (11) 

B B 

llGir ,7,, ,_ , 227r //t,x2 /- \2n 

"'' SIM^ W(^'S) + 27jt/4 OLeedes{{dd) + (ss) ), 

/i±e'^ = -^2(1 + — + e«e,^)M|(l-e"%)--^(m„+m,)2lnAfs(l-e"^) 

47r TT TT 47r 

-^("^n(ss)+m,(Su)) + ^2lf|<V^'^ (12) 

1767r \ /- \ 227r o /- \2n 

+ gY]g4"«W(««) + 27M4"ee„e,((uu) + (ss) ), 

fhoe'^ = -^(l + ^ + e,e,^)M|(l-e"^)--^(md-m,)2lnMB(l-e"^) 

+ ^(-.(^-^> + jin^i^G') - ^aAdd)i,s) (13) 

B B B 

B B 

fh±e^^ = -^{l + — + e^es — )Ml{l~e^^)~^{mu-msf\nMB{l-e~^) 

47r TT TT 47r 

+ ^Ti^u{ss) + ms{uu)) + J^^i-^C!^) - ^^o:s{uu){ss) (14) 



III. NUMERICAL RESULTS 

For numerical calculations, we use the following values of condensates [TIL Il3. ITHl ITfil ItR ITM Il9l | : 

= ^^^^^^ ^ -(0.240 GeV)^ 

A = - 1 = -0.0074 

[uu) 

(ss) = (0.8 ±0.1) X (15) 
{—G^) = 0.012 GeV'^, 

TT 

m„ = 5.3MeV, ma = 9.4MeV, = ISOMeV, 

ap — , a, = 0.7. 

137 

The up and down quark condensates have uncertainly in the absolute values. However, we keep their difference 
A w -0.0074. 



A. The QCD Sum Rule for the K meson 

Differentiating Eqs. (|12|) and (|13|l with respect to and dividing the results by themselves, we obtain the masses 
for K° {K°) and mesons. For the study of the QCD sum rule, we have two parameters, the threshold value sq 
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and the Borel mass Mb- Herein below we study the Borel mass dependence in the region 0.5 < Mg < 3.0 GeV^ with 
So ^ 0.9 GeV^. Further discussions on the Sq dependence wiU be presented in the end of this work. 

For the absolute values of the mass of the K meson, the present QCD sum rule does not work well, because K is 
the Nambu-Goldstone boson having a strong collective nature due to the non-perturbative QCD dynamics. However, 
the isospin symmetry breaking effects can reasonably be studied in the QCD sum rule. 

The mass difference {Auik — mxo(K°) ~ ^k±) is shown in Fig. [3 as a function of the Borel mass square Mg. 
The dashed curve is obtained when the threshold value sq = 0.900 GeV^ is used both for {K^) and K^. The 
resulting mass difference turns out to be negative which does not agree with the experiment. Also the Borel stability 
is not good. We can fine tunc the threshold value sq and use different values for (K^) and K^. The solid line 
is obtained when we take so{K°,K°) = 0.916 GeV^ and sa{K^) = 0.900 GeV^, with which the sum rule value takes 
ArriK = 4± 1 MeV (Af| > 1 GeV^). This is consistent with the experimental value Atuk = 3.972 ±0.027 MeV [T^. 
The Borel stability is also improved for Mg > 1 GeV^. 




Borel Mass^ [GeV^] 

FIG. 2: The mass difference of the K meson, as a function of the Borel mass square M|. The dashed curve is obtained when 
so{K",K°) = so(K^) = 0.900 GeV^ The solid curve is obtained when so{K" , K°) = 0.916 GeV^ and so{K^) = 0.900 GeV^ 



Now let us study the K decay constant. We need to in put the mass of the K meson which we use the experimental 
values, TO^o(;fo) = 497.6 MeV and mj^± = 493.7 MeV The left panel of Fig. |3 shows the K"^ decay constant 
/if± as a function of the Borel mass square when Sq = 0.900 GeV^ is used. The result for K'^ {K^) and can 
not be distinguished in this figure (see the right panel and discussion below) . It is interesting that the sum rule values 
take around 165 MeV with a good Borel stability and is consistent with the experimental value fi(± — 159.8 ± 1.84 
MeV 14]. 

The difference of the K decay constants A/k = fK"{KO) ~ /-ftr± is plotted in the right panel of Fig. O as a function 
of the Borel mass square M^. The meaning of the dashed and solid curves are the same as for Fig. |21 When the 
same threshold values are used, AJk takes values 0.2 ~ 0.5 MeV for 1 < < 3 GeV^ with some strong Borel mass 
dependence. However, when using the different threshold values, we obtain A/x — 1.5 ± 0.2 MeV with a good Borel 
stability for 1 < Af| < 3 GeV^. 



B. The QCD Sum Rule for the K* Meson 

For the K* meson, we expect that the QCD sum rule works well just as in the case of the p meson. In order to 
check the validity of the present sum rule, we show the mass of K*^ in the left panel of Fig. ^ where we find a very 
good Borel stability. The absolute value depends slightly on the choice of the threshold value sq, which we choose 
So = 1.80 GeV^ to reproduce the experimental value m^^.i = 891.7 MeV jQ. The result for {K*°) is very 
similar. 

The mass difference ArriK- = to^*o — mj^'.± is shown in the right panel of Fig. ^ as a function of the Borel mass 
square Mg. The dashed curve is obtained when the same threshold value sq = 1.80 GeV^ is used both for K*° (K*^) 
and K*^. The Borel stability is not good. We can fine tune the threshold value so again and use different ones for 
K*° {K*°) and K*"^. The solid line is obtained when we take so{K*°,K*°) = 1.83 GeV^ and soiK*"^) = 1.80 GeV^, 
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3: The left panel shows the decay constant fx±, as a function of the Borel mass square 
So — 0.900 GeV^. The right panel shows the difference of the K decay constants, as a function of the Borel mass square Mg 
dashed curve is obtained when so{K°, K°) = so{K^) = 0.900 GeV^ 
- 0.900 GeV^ 



FIG. 



and so{K 



for threshold value 
The 

The solid curve is obtained when so{K'',K°) = 0.916 GeV^ 



with which the sum rule value takes ArriK* = 7 ± 1 MeV {Mg > 1 GeV^). This is consistent with the experimental 
value AniK* = 6.7 ± 1.2 MeV The Borel stability is much improved for Af| > 1 GeV^ 
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FIG. 4: The masses of K*°{K*°) and K*"^ and their differences, as a function of the Borel mass square Mg for threshold values 
s^.o R'o = 1-83 GeV^ and Sj^,± = 1.80 GeV^ 



For the study of K* decay constant, we use the experimental values for the mass of K* meson, mj^to = 896.1 MeV 
and mx*± — 891.7 MeV [lj|. The left panel of Fig. Elshows the K*^ decay constant /j<-.± as a function of the Borel 
mass square Mg when sq = 1.80 GeV^ is used. The result for (K*^) is very similar. The sum rule values take 
around 230 MeV with a good Borel stability. We can estimate the decay rate of r~ K*'Vr as j2Qj 

= 3.42 X 10"" GeV, 

which is not far from the experimental values Tt-^K'-u^ = 2.78 x lO"^"* GeV [T3 |. 
We can also calculate the decay rate oi ^ Vr 

■ = |''"-|"S(7£r)""«('-t)' 

= 1.67 X 10"" GeV, 

which is also consistent with the experimental values T^-^i^-^^ — 1.48 x 10"^'' GeV [l^. 
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The difference of the K* decay constants Afx* = fK'°(K'°) ~ Jk'^ is plotted in the right panel of Fig. |S1 as a 
function of the Borel mass square . The meaning of the dashed and solid curves are the same as for Fig. ^ 




0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 

Borel Mass^ [GeV^ ] Borel Mass^ [GeV^ ] 

FIG. 5: The decay constants of K*'^{K*'^) and K*^ and their differences, as a function of the Borel mass square Mg for 
threshold values Sjf,ajf,a — 1.83 GeV^ and s^*± — 1.80 GeV^. 



C. The Threshold Value so Dependence 

Finally, let us investigate sq dependence of the present QCD sum rule analyses. In order to see its typical behavior, 
we fix the Borel mass to be Mg — 1 GeV^. We use the two threshold values 

so{K", K*°) = So + Aso, for K" or K*° , 
sq{K^,K*^) = So, for iC^ or is:*± . 

As So is varied, Asq is fixed such that the experimental mass difference Attik or ArriK* is reproduced. The differences 
of the decay constants A/ are then computed as functions of sq. The resulting Aso and A/ are plotted in Fig. for 
K and in Fig. [7| for K*. It is interesting to observe that although Aso are monotonically increasing functions, A/'s 
are rather stable as sq is varied. It would be an indication that the present sum rule analyses especially for A/ are 
stable. 




0.7 1 1.5 2 2.5 3 0.7 1 1.5 2 2.5 3 

So [GeV^] So [GeV^] 



FIG. 6: Aso = so(A' ) — so (A" ) and A/k as functions of sq. Aso is determined so as to reproduce the experimental value of 
ArriK ^ 3.972 MeV at M| = 1 GeV^ 



IV. SUMMARY 



In this paper, we have studied isospin breaking for masses and decay constants of K and K*. We have adopted 
gauge invariant currents coupled by a photon field. We have then estimated isospin symmetry breaking effects 
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FIG. 7: Aso = So (if*") — so(-?C'* ) and A/k* as functions of sq. Aso is determined so as to reproduce the experimental value 
of Atuk- = 6.7 MeV at M| = 1 GeV^ 



through different values of the parameters such as quark masses, condensates and threshold values. Quark masses 
and condensates were fixed from other sources, while the threshold values were fixed such that the mass differences of 
charged and neutral K and K* were reproduced. The resulting decay constants were found to be very stable against 
the change in the Borel mass and the threshold values. The resulting values for Am and A/ are consistent with 
experimental values. 

The present analysis with good stability indicates that the QCD sum rule can be applied to study the symmetry 
breaking effects in hadron physics. In the near future, BESIII collaboration will measure the mass splittings of K and 
K* systems precisely. Investigation of isospin symmetry breaking patterns helps to explore the low-energy sector of 
the underlying QCD dynamics. 



Acknowledgments 



H.X.C is grateful to the Monkasho fellowship for supporting his stay at Research Center for Nuclear Physics 
where this work is done. A.H. is supported in part by the Grant for Scientific Research ((C) No. 16540252) from 
the Ministry of Education, Culture, Science and Technology, Japan. S.L.Z. was supported by the National Natural 
Science Foundation of China under Grants 10375003 and 10421503, Ministry of Education of China, FANEDD, Key 
Grant Project of Chinese Ministry of Education (NO 305001) and SRF for ROCS, SEM. 



T. Das, G. S. Guralnik, V. S. Mathur, F. E. Low and J. E. Young, Phys. Rev. Lett. 18, 759 (1967). 

K. Fujikawa and R J. O'Donnell, Phys. Rev. D 9, 461 (1974). 

L. Scorzato, Eur. Phys. J. C 37, 445 (2004) arXiv:hep-lat/0407023 . 

W. A. Bardeen, J. Bijnens and J. M. Gerard, Phys. Rev. Lett. 62, 1343 (1989). 

A. Duncan, E. Eichten and H. Thacker, Phys. Rev. Lett. 76, 3894 (1996) arXivih lplat /9602005| . 
S. L. Zhu and Z. P. Li, Phys. Rev. D 55, 7093 (1997) arXiv:hep-ph/9703371 . 

L. S. Kisslinger and Z. P. Li, Phys. Rev. Lett. 74, 2168 (1995) arXiv:hep-ph/9409377 . 

M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Nucl. Phys. B 147, 385 (1979). 

P. Ball and R. Zwicky, Phys. Lett. B 633, 289 (2006) arXiv:hep-ph/0510338 . 

L. J. Reinders, H. Rubinstein and S. Yazaki, Phys. Rept. 127, 1 (1985). 

J. Gasser and H. Leutwyler, Nucl. Phys. B 250, 465 (1985). 

T. Hatsuda, H. Hogaasen and M. Prakash, Phys. Rev. C 42, 2212 (1990). 

V. L. Chernyak and A. R. Zhitnitsky, Phys. Rept. 112, 173 (1984). 

S. Eidelman et al. [Particle Data Group], Phys. Lett. B 592 (2004) 1. 

K. C. Yang, W. Y. P. Hwang, E. M. Henley and L. S. Kisslinger, Phys. Rev. D 47, 3001 (1993). 

V. Gimenez, V. Lubicz, F. Mescia, V. Porretti and J. Reyes, Eur. Phys. J. C 41, 535 (2005) arXiv:hep-lat/0503001| 
M. Jamin, Phys. Lett. B 538, 71 (2002) arXiv:hep-ph/0201174j^ 

B. L. lofle and K. N. Zyablyuk, Eur. Phys. J. C 27, 229 (2003) |arXiv:hep-ph/0207183| . 



10 



[19] A. A. Ovchinnikov and A. A. Pivovarov, Sov. J. Nucl. Phys. 48, 721 (1988) [Yad. Fiz. 48, 1135 (1988)]. 

[20] J. F. Donoghue, E. Golowich and B. R. Holstein, Dynamics of the Standard Model (Cambridge University Press, 1992) 



